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ABSTRACT 

We construct high signal-to-noise "template" spectra by co-adding hundreds of spectra of nearby 
dwarfs spanning K7 to M4, taken with Keck/HIRES as part of the CaHfornia Planet Search. We 
identify several spectral regions in the visible (370 - 800 nm) that are sensitive to the stellar luminosity 
and metallicity. We use these regions to develop a spectral calibration method to measure the mass, 
mctallicity, and distance of low-mass stars, without the requirement of geometric parallaxes. Testing 
our method on a sample of nearby M dwarfs we show that we can reproduce stellar masses to about 
8 - 10%, metallicity to ^ 0.15 dex and distance to 11%. We were able to make use of HIRES spectra 
obtained as part of the radial velocity monitoring of the star KOT314 to derive a new mass estimate 
of 0.57 ±0.05 Mq, a radius of 0.54 ±0.05 i?0, a metallicity, [Fe/H], of -0.28 ±0.10 and a distance of 
66.5 ± 7.3 pc. Using HARPS archival data and combining our spectral method with constraints from 
transit observations, we are also able to derive the stellar properties of GJ 3470, a transiting planet 
hosting M dwarf. We estimate a mass of 0.53 ± 0.05 Mq, a radius of 0.50 ± 0.05 Rq, a metallicity, 
[Fc/H], of 0.12 ± 0.12 and a distance of 29.9±j^:| pc. 



1. INTRODUCTION 

M-type dwarf stars are poorly understood compared 
to higher-mass FGK stars because of the difficulty in 
modeling both their atmospheres and interior stru ctures 
(|Hauschildt et allHOOl IChabrier fc Barafrdl2000l \. The 
cool atmospheres of these stars contain many molecu- 
lar species such as VO and TiO, which dominate both 
the line and continuun i opacity in their photospheres 
(jKirkpatrick et aLlll991[ ). Due to the millions of molec- 
ular transitions, many of the opacity sources have yet 
to be accounted for, making it difficult to synthesize M- 
dwarf spectra. Modeling has also proven challenging in 
this mass regime because convection plays an important 
role in the structure of the star ([Baraffc ct al. 1998; Ribas 
120061 [Morales et aIll2010D . Much of the input physics for 
low-mass stars remains approximate, and consequently, 
the physical properties of M dwarfs outside of binary sys- 
tems are difiicult to measure (for example studies of bi- 
narie s see LoDCz-Moralcs & Ribas 2005; Coughli n et al.l 

[mil). 

Despite these challenges, there have been many at- 
tempts to discern the intrinsic properties of isolated M 
dwarfs. These efforts have focused on using observational 
properties to either aid the theoretical treatments or to 
develop empirical calibrations based on stars of known 
physical properties. iDelfosse et a l. (2000) demonstrated 
a tight relation between absolute magnitudes in the in- 
frared passbands (J, H, and Kg) and stellar mass as 
determined from binary systems. As the metallicity of 
these stars increases, their colors redden and the bolo- 
metric luminosity decreases. These effects balance in the 
infrared to produce relatively tight mass-luminosity re- 
la tions. 

IDelfosse et al.l ()2000f ) also suggested that the large 

^ Corresponding author: jspineda@astro.caltech.edu 
^ California Institute of Technology, Department of Astron- 
omy, 1200 E. Cahfornia Ave, Pasadena CA, 91125, USA 



scatter in the color- magnitude d iagram, (V — Kg) - 
Mks, was due to this same effect. iBonfils et all (|2005aD 
took advantage of this sensitivity to metallicity to de- 
velop a photometric calibration based on wide bina- 
ries with an M-dwarf secondary an d an FG K primary 
of known [Fe/H]. The IBonfils et al.l ()2005aD photomet- 
ric calibration has since been revised using improved 
V-ba nd photometry ([Johnson k, Appsll2009l iNeves et al.l 
1 201-^) and by using a phys i cal model based guide 
( Sch Jaufman fc Laughlin 2010*). Schl aufman fc Laughlinl 
(i201Ql) and Ncvcs ct al. (2012) parameterize the metallic- 
ity as a function of /S.{V — Ks)^ which measures how much 
a given star deviates from the main sequence. These ef- 
forts represent significant progress in determining the in- 
trinsic properties of M dwarfs; however they require pre- 
cise photometry and a parallax measurement. Although 
these requirements are met for nearby M dwarfs (d < 15 
pc), the lack of parallaxes for fainter, more distant stars 
severely limits our knowledge of M dwarfs beyond the 
Solar neighborhood. 

Additionally, spectral synt hesis has been used to mea- 
sure M-dwarf me t alhcities ([Woolf fc WallersteinI 120051 : 
iBean et al.l 120061: lOnehag et all I2012D . For higher- 
mass stars, spectral line comparisons based on equiv- 
alent width measurements can be used to compute 
abundances. However, the line blanketing in low-mass 
stars makes this extremely difficult, as individual lines 
bl end and cornplete ly dominate any thermal continuum 
. IBean et all ([2006D applied atmospheric models to bi- 
nary systems composed of an FGK primary and an M- 
dwarf secondary. Assuming a co-evolutionary system, 
the metallicity of the primary can be determined accu- 
rately and used t o calibrate the mo deling of M-dwarf 
spectra. However, IBean et aD ([2006f ) were off by ~ 0.3 
dex compared to later photometric calibrations, likely 
due t o deficiencies in the low-ni ass atmospheric mod- 
els. IWoolf fc WallersteinI ([20051 ) were successful, but 
only because they used extremely metal-poor M dwarfs 
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with minimal molecular line blanketing. More recently, 
lOnehag et al.l (j2012D matched synthetic spectra to their 
high-resolution J-band spectra to measure metallicities 
to within 0.09 dex of the best photometric calibrations. 

Instead of relying on atmospheric models, recent, more 
accurate met hods have made u s e of o bservation-based 
calibrations. iRoias-Avala et all (|2010( ) used moderate- 
resolution infrared spectra to develop metallicity indica- 
tors b ased on Ca I and Na I features in the i^-band. Sim- 
ilarly, iTerrien et al.l (|2012f ) used Ca I and K I features in 
the iJ-band to measure metallicity. These near-IR spec- 
troscopic calibra tions agree well with current photomet- 
ric c alibrations ()Roias-Avala et all I2012t ITerrien et al.l 
[20TI . 

Accurately measuring the properties of low-mass stars 
has gained renewed urgency because of the discovery of 
a multitude of planets around M dwarfs (iButler et al.l 
[200llRivera et al.ll2005l: iBonfils et al.ll2005bD . There are 
also many M dwarfs among the hosts of planet candi- 
dates discovered by the Kepler Mission fBoruck i et al.l 
[20Tlt IBatalha et al.|[20T2l : IMuirhead et a l. 2012a). Ad- 
ditionall y, these stars may host a multitude of terrestrial 
planets (IMuirhead et all l2012bl: iBuchhave eFall 120121: 
iSwift et al.l |2013() " For radial velocity and transit de- 
tected planets, determining the physical properties of the 
planets requires an accurate measurement of the physical 
properties of their host stars. 

In this contribution we develop a method to measure 
the physical properties of low-mass stars using spec- 
troscopic indices from high-resolution optical spectra. 
Specifically, our technique provides estimates of the ab- 
solute NIR magnitudes (TWj, Mh, Mk), distances d, 
and A(y — Kg) for M dwarfs without parallaxes. Us- 
ing known observational calibrations these quantities 
can be converted to est i mate both mass and metallic- 
itv (iDelfosse et all I2OOOI: iSchlaufman fc LaughlinI I20T0I: 
INeves et al.ll2012D . We create a library of high signal-to- 
noise, high-resolution "template" low-mass dwarf spectra 
with known photometric properties to develop calibra- 
tion curves based on the strength of various absorption 
features. By measuring the strength of these features 
in the spectrum of a star with unknown properties and 
using the calibration curves, the photometric properties 
of the low-mass star can be estimated which in turn can 
be used to estimate the star's physical properties. In 
Section 2, we describe the data sample used to develop 
our library. In Section 3, we cover our analysis meth- 
ods to construct and develop our calibration curves. In 
Section 4, we apply our methods to a distant M-dwarf 
listed among the host stars containing Kepler exoplanet 
candidates and another star with a recently discovered 
transiting neptune-mass planet. Lastly, in Section 5, we 
discuss the utility of our methods. 

2. DATA 
2.1. Sample 

Over the past decade the California Planet Sur- 
vey (CPS) has obtained spectroscopic measurements of 
more than 2,500 stars at Keck Observatory, monitor- 
ing their radial velocities for the characteristic signa- 
ture o n the host star ind uced by the presence of a 
planet (jHoward et al.|[20Tol) . We make use of thei r High- 
Resolution Echelle Spectrometer (HIRES ; iVogt et al.l 
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Fig. 1. — Distribution of stellar properties in sample. Top panel 
is A^/fs, and bottom is A(V — Ks)- The calibration sample thus 
spans a broad range or stellar properties 

119941 ) observations of 155 M dwarfs to build a high- 
resolution, high signal-to-noise "template" spectrum of 
each star. The typical CPS program HIRES set-up gives 
a resolving power R = A/AA « 50,000 and uses an 
iodine-cell for precise wavelength calibrati on of radial ve- 
locity measurements (jHoward et al.l[20To[ ). In our appli- 
cation the iodine lines are contaminants, so we can only 
make use of the blue and red chips of the HIRES detec- 
tor, where there are no lines from the iodine cell. The 
blue portion of the spectrum allows us to examine be- 
tween 370 nm and 480 nm, while the red portion of the 
spectrum covers between 650 nm and 800 nm. 

The M dwarfs in the CPS sample have well de- 
fined 2MASS p hotometry and parallaxes from Hippar- 
cos (jPerrvman fc others. 1997; Cutri et al. 2003). This 
allows us to characterize the stellar sample in terms of 
their absolute near-IR magnitudes, {Mjhks) and color, 
{V ~ Kg). We culled stars from the sample that were 
overly active or young, selected on the basis of hav- 
ing published rotation periods less than 5 days, high 
X-ray luminositi es with a X-ray co unt rate > 1 count 
in ROSAT ()Voges et al.l 1 19991) or as being mem- 
bers of a young open cluster or moving group. We also 
culled stars that were known to be unresolved binaries 
or turned out to be missing geometric parallax measure- 
ments. We also limited the sample to stars brighter than 
■Mks = 8, so that the range of properties is well sam- 
pled by the CPS stars. This left us with 119 calibra- 
tion stars. The sample spans a range of A4ks from 
4.5 to 7.5 as shown in t he to p panel of Figure [TJ Us- 
ing the lJohnson fc A pps' (2009") solar-metallicity main se- 
quence relation, defined as Mks ~ J2 ai(^~^s)% where 
a = {-9.58933,17.3952,-8.88365,2.22598,-0.258854, 
0.0113399}, we can calculate the quantity, A{V — Kg), 
the difference between the observed color and the main- 
sequence color for a star of the same A4ks such that 
positive values of A{V — Kg) correspond to redder ob- 
jects. We will henceforth refer to /S.{V — Kg) as the 
color offset, which can be used as a proxy for metal- 
hcity (jSchlaufman fc Laughlinll2010D . The bottom panel 
of Figure [1] shows that our sample spans a broad range 
in color offset and hence a broad range in metallicity. 
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2.2. Spectral Library 

Since our sample of calibration stars have been moni- 
tored for radial velocity shifts indicative of planets over 
the past 10-15 years, each star has an average of 25 ob- 
servations. By combining the individual observations we 
can produce high signal-to-noise, high-resolution tem- 
plate spectra for each of the calibration stars in our sam- 
ple. We rebinned each star's spectrum onto a wavelength 
scale that is evenly spaced in In A (cadence of 1.9 x 10^^), 
which allows us to properly Doppler shift the sp ectra 
with respect to one another ()Tonrv fc Davis! 119791) . We 
also corrected for small differences in the wavelength so- 
lution from night to night (on the order of a couple of 
pixels) for every spectral order due to changes between 
the cross-disperser angle and the echelle angle as well as 
changes in the slit illumination for any given observation. 
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Fig. 2. — Normalized Intensity as a function of wavelength for 
sensitive regions with bin-averaged spectra of width 0.1 mags in 
A4ks- Colors correspond to Mks in the range 4.5 to 7.5, with red 
being fainter stars and blue being brighter. 

Having aligned the spectra and removed defects (cos- 
mic ray hits etc.), we simply co-added the flux to produce 
a high signal-to-noise ratio spectrum for each star. The 
red portions of the spectrum yield a typical signal-to- 
noise ratio (SNR) per resolution element ranging from 
~ 800 to ~ 4000, depending on the number of observa- 
tions for the particular star. For the blue-chip spectra 
the SNR ranges from - 100 up to - 600. The SNR 
of the red side is higher primarily because the peak of 
the M dwarfs' spectral energy distribution lies closer to 
near-infrared wavelengths than to the blue portion of the 
spectrum. 

In order to compare the spectra it was necessary to 
normalize each order to remove the effect of the blaze 
function of the spectrometer and to account for a pseudo- 
continuum. We normalize each spectral order individu- 
ally. First, we separate the order into tens of "chunks", 
masking out telluric regions, and take the top 1-2% of 
points of each chunk as representing the "continuum" . 
We then fit a low-order polynomial to the continuum 
points across the full extent of the order. Lastly, we 
divide the spectral order by this polynomial, ignoring 
problematic points at the ends of the orders, to get the 
normalized spectrum (see Appendix for further details). 
Since the echelle spectra are not flux calibrated and the 
blaze function distorts the shape of the spectra, a pseudo- 



continuum is not well defined through our normalization 
procedure. However, because we use the same contin- 
uum regions for all stars, this process allows for reliable, 
differential comparisons among stars of different types. 

In Figure [U we plot the template spectra averaged in 
A4ks bins with width of 0.1 mag. The different colors 
correspond to different values of Mks with red corre- 
sponding to cooler stars and blue to hotter stars. In 
the Figure, we see that the absorption features are quite 
distinct from the "continuum" , which match across the 
different stars. It is clear in the high signal-to-noise spec- 
tra that features such as these are quite sensitive to Mks 
and we can use the strength of the absorption for a given 
spectrum as indicative of the stars intrinsic luminosity. 
Deviations from a strict sequence are primarily due to 
metallicity effects. At a given Mks changes in metallic- 
ity will affect the strength of certain absorption features. 
Accounting for this second-order effect will provide us 
with a valuable indicator of stellar metallicity. 

3. ANALYSIS 
3.1. Spectral Calibration 

The presence of spectral regions sensitive to changes 
in MjHKs motivates the development of a quantitative 
relationship between the strength of each region and the 
physical properties of a star. The regions that were found 
to be sensitive and useful for calibration are listed in Ta- 
ble [T] We identified these regions by eye looking across 
the full spectrum, ignoring telluric regions, and requir- 
ing the continuous regions of the spectrum to have have 
monotonically increasing absorption with decreasing stel- 
lar effective temperatures (see Figure [2]). The useful re- 
gions consist predominately of portions of TiO and VO 
bands. Note that the spectral indices listed in the ta- 
ble are all in the red portion of the spectrum. Although 
there also appeared to be sensitive regions at blue wave- 
lengths, the lower overall flux level limits their usefulness, 
both for the calibration procedure and for future obser- 
vations. In addition to the regions listed in Table [T] there 
were other regions that we identified as sensitive to the 
physical stellar properties. However, we selected a sub- 
set that when combined provided the optimal calibration 
(see Section [231) • 

The CPS sample spans a representative range of prop- 
erties for low-mass stars, making it useful for our cali- 
bration procedure. Since the sample is limited in size 
and individual stars only represent discrete points in the 
mass-metallicity plane, comparing the spectra directly 
is less than ideal, giving poor parameter resolution. In- 
stead, it is preferable to compare the strength of sensitive 
features, measured from their integrated fluxes (equiva- 
lent widths, EWs), allowing us to flt smooth functions 
to observed trends in EW and providing a continuous 
relationship between the integrated flux and the stellar 
properties. The integrated flux is defined as 

EW = AX- J S{X)dX (1) 
^AX-Y,S^hx, (2) 

i 

where AA is the width over which the integral is com- 
puted and S'(A) is the normalized spectrum as a function 
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TABLE 1 
Spectral Indices 



Center 


Total 


Molecule / 


Wavelength (nm) 


Width (nm) 


Line 


658.24 


0.820 


TiO 


uuu. oo 


0.59 


TiO 




0.81 


TiO 


uuo. 


1.71 


TiO 


710.27 


2.91 


TiO 


713.64 


2.48 


TiO ^ 


727.01 


0.06 


TiO 


770.04 


0.79 


TiO / KI 


770.83 


0.50 


TiO 


776.92 


2.38 


TiO 


780.13 


0.92 


TiO 


787.30 


0.85 


TiO / VO 


791.04 


0.51 


VO 


792.26 


1.67 


VO 


793.26 


0.30 


VO 


794.87 


0.24 


VO 



Predom inant Molecules and L ines based on 
spectra of IKirkpatrick et al.lll99ll 
^ This region is closely related to the Ti04/5 
bands defined in lHeid et al.|[T995l 

of wavelength A. Equation [2] gives the approximation 
for discretely sampled spectra over pixels that span AA 
evenly sampled with width h\. To determine the er- 
rors on our equivalent width measurements we randomly 
simulate the spectral observation using poisson statis- 
tics with a mean in each spectral bin given by the pho- 
ton counts of the actual data. We take the error on the 
equivalent width measurement to be the standard devi- 
ation of the distribution of EW values in the simulation. 
The typical error on the EW measurements are ~ 1 — 5%. 




(V-K) 



Fig. 3. — EW plotted as a function of A4ks for all stars in the 
sample. The colors are ordered according to the A{V — Ks) of the 
stars in the sample. There is a clear gradient that corresponds to 
differences in metallicity. The lines correspond to contours in our 
polynomial model of constant color offset, A{V — Ks), with values 
0.5, 0.0 and -0.5 going from top to bottom respectively. This region 
is centered at 770.04 nm with a width of 0.50 nm 

In Figure[3]we plot the behavior of EW as a function of 
Mks and A{V—Ks) for a particular spectral region. The 
colors correspond to A{V — Kg) with red corresponding 
to stars that are more metal-rich and blue corresponding 
to stars that are more metal-poor. For spectral regions 
such as this one, the strength of the feature increases with 
increasing metal content as well as decreasing luminos- 
ity. The behavior, expressed in Figure [3l motivates the 



parameterization of each EW in terms of M.jhks and 
/S.{V — Kg). For a given region, Z, we calibrate the equiv- 
alent widths against each of the absolute magnitudes and 
the color offset. Our calibration EW, which allows us to 
interpolate between the discrete sample star properties, 
is given by 



EWz,, 



2 

E 

i=0 



(3) 



where a G { J, H, Ks }, so there is a separate calibra- 
tion for each passband using the same spectral region. 

We fit each passband separately instead of going di- 
rectly to mass, despite tight mass-luminosity relations, 
because the infrared colors (ex. J — A') are not simply 
functions of mass and can change with metallicity; this 
is in contrast to other broadband ph o tometr ic studies of 
low- mass stars (j Johnson et al.l l20lTl I2012bl ). In addi- 
tion to the simple polynomial terms expressed in Equa- 
tion [21 we include a cross term governed by the coef- 
ficient c that accounts for differences between brighter 
and fainter stars in how their absorption strength re- 
sponds to the addition of metals. The need for such a 
term is evident in how TiO features are known to satu- 
rate in late-type M dwarfs and that VO features are not 
appare nt in early M dwarfs bu t appear in late-type M 
dwarfs (jKirkpatrick et al.lll991[ ). 

3.2. Fitting Broadband Photometry 

We used a Bayesian method to fit Equation |3] to each 
set of EWs. In addition to the coefficients in Equation|3l 
we also incorporated an additional parameter, ct, to take 
into account intrinsic scatter in our choice of parameter- 
ization. Using a Markov Chain Monte Carlo (MCMC) 
technique with a Metropolis-Hastings algorithm, we ex- 
plore the posterior distribution for the best parameters 
in the calibration conditioned on the known properties 
of the stars in the sample. The best calibration param- 
eters are those that maximize their respective marginal 
distributions, and thus maximize the probability of each 
parameter reproducing the stellar properties of the stars 
in the sample with our simple model. We report the pa- 
rameter values of our calibration in Table [2] For a given 
spectral index the first row of the table corresponds to 
the parameters for the calibration with M.j, the second 
row for Mh and the third for Mks- 

These curves provide a calibration, which we can apply 
to stars of unknown properties with measured EWs and 
estimate AAjhks and A(V^ — Kg). Taking the errors in 
the integrated fiux as normally distributed, we minimize 
a statistic to get the stellar photometric properties. 
The statistic is split up into a couple components. 



X 



E2 2 
Xi + Xp , 



(4) 



where the first term comes from the spectral calibration 
and is summed over the different indices, /, and the last 
term correspond to photometric constraints. The sum 
over all the different indices allows us to partially break 
the degeneracy between metallicity and mass inherent in 
the individual indices. 
For each spectral region, I, the fitting statistic is 
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TABLE 2 

Calibration Parameters 



Center (nm) 


10 (nm) 


ai (nm) 


02 (nm) 


b (nm) 


c (nm) 


cr (nm) 


658.24 


-0.4501058 
-0.3490226 
-0.3486802 


0.09198486 
0.07719585 
0.07969309 


-0.001427114 
-0.0005361921 
-0.0005093491 


0.07882096 
0.06294231 
0.07282996 


0.002933661 
0.006100980 
0.005513314 


0.01270252 
0.01208604 
0.01123767 


660.35 


-0.1596480 
-0.1520269 
-0.1553699 


0.01916057 
0.02420063 
0.02603111 


0.001997996 
0.001671400 
0.001767160 


0.001896713 
0.003454885 
0.01073864 


0.008210727 
0.009034893 
0.008722165 


0.007939490 
0.007594586 
0.007079715 


663.89 


-0.4226041 
-0.3302306 
-0.3803160 


0.08017676 
0.06765508 
0.08669426 


-0.00008847487 
0.0007076964 
-0.0005393973 


0.07432834 
0.05919351 
0.08166828 


0.003435445 
0.006587219 
0.004118372 


0.01253907 
0.01168440 
0.01082510 


668.82 


-2.836934 
-2.196094 
-2.374936 


0.7375158 
0.6257768 
0.7089150 


-0.03934088 
-0.03455122 
-0.04205847 


0.6871808 
0.5937715 
0.6733981 


-0.06801754 
-0.05868804 
-0.07173435 


0.03053437 
0.02922317 
0.02711429 


710.27 


-5.145841 
-4.044212 
-4.307040 


1.292010 
1.116565 
1.244481 


-0.06165890 
-0.05456582 
-0.06565496 


1.147429 
1.010953 
1.162449 


-0.08640261 
-0.07058973 
-0.09311087 


0.06357273 
0.06100981 
0.05549368 


713.64 


-6.357848 
-5.054271 
-5.253722 


1.673048 
1.471289 
1.593891 


-0.08776160 
-0.08102439 
-0.09263240 


1.399011 
1.272483 
1.423004 


-0.1214399 
-0.1101706 
-0.1339820 


0.05989609 
0.05812801 
0.05169340 


727.01 


-0.1138363 
-0.09486709 
-0.09856623 


0.02944591 
0.02712001 
0.02938447 


-0.001551650 
-0.001531907 
-0.001748182 


0.02305361 
0.02163241 
0.02431960 


-0.001910082 
-0.001814339 
-0.002235494 


0.001294007 
0.001242313 
0.001166170 


770.04 


-0.1598857 
-0.1181027 
-0.1085764 


0.008933704 
0.008006942 
0.004873543 


0.006038753 
0.006416439 
0.007252633 


0.1399386 
0.1134480 
0.1158942 


-0.004808356 
-0.0005393014 
-0.000001015819 


0.01121777 
0.01118438 
0.009809055 


770.83 


-0.3815610 
-0.2753433 
-0.3191640 


0.07699559 
0.05733777 
0.07383029 


-0.001319326 
-0.00003557394 
-0.001163489 


0.04795545 
0.03629728 
0.05463092 


0.003966268 
0.006487797 
0.004446250 


0.009098964 
0.008746648 
0.007962008 


776.92 


-1.408374 
-1.204047 
-1.302366 


0.2440679 
0.2346570 
0.2751659 


0.001240992 
0.001405473 
-0.0008869298 


0.1658279 
0.1440118 
0.2068376 


0.02825762 
0.03585729 
0.03001058 


0.04063114 
0.03961519 
0.03584416 


780.13 


-0.3511850 
-0.2988942 
-0.2972508 


0.05775816 
0.05584809 
0.05705502 


0.001359554 
0.001437870 
0.001681733 


0.05715288 
0.05035518 
0.06221154 


0.007504470 
0.009766074 
0.009053781 


0.01079916 
0.01034952 
0.009174584 


787.30 


-0.07969772 
-0.04158101 
-0.06175869 


-0.01437250 
-0.02027625 
-0.01473584 


0.005574401 
0.006402968 
0.006440116 


0.009402555 
-0.003833748 
0.005074938 


0.01204589 
0.01558149 
0.01547640 


0.01161012 
0.01126873 
0.01026771 


791.04 


0.01844099 
0.01743679 
-0.0009558277 


-0.02735590 
-0.02505636 
-0.02019831 


0.004700167 
0.004849061 
0.004791414 


-0.02789386 
-0.02857221 
-0.02164277 


0.01186093 
0.01330172 
0.01312965 


0.006290352 
0.006062947 
0.005528097 


792.26 


0.2785090 
0.2310036 
0.2100510 


-0.1369218 
-0.1224368 
-0.1211605 


0.01704969 
0.01723172 
0.01818554 


-0.1025768 
-0.09838984 
-0.09387771 


0.03602162 
0.03933845 
0.04123889 


0.01778133 
0.01729166 
0.01595478 


793.26 


0.1454432 
0.1161707 
0.1180228 


-0.04912226 
-0.04273733 
-0.04515902 


0.004548514 
0.004429765 
0.004855676 


-0.01954720 
-0.02024410 
-0.02172495 


0.006100179 
0.006868234 
0.007529807 


0.003336126 
0.003257263 
0.003095342 


794.87 


0.04278556 
0.03241636 
0.02124669 


-0.02045305 
-0.01713182 
-0.01420778 


0.002583565 
0.002513652 
0.002431369 


-0.004188786 
-0.002884554 
-0.0002786377 


0.004094394 
0.004370116 
0.004293440 


0.002292025 
0.002306728 
0.002088636 



^ The first row for a given index center corresponds to parameters for the J cahbration, the second row for H 
caUbration and third row for Ks calibration 
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can use the distance, d, as an additional parameter by 
requiring that our estimates of Mjhks reproduce the 
observed photometry. 



where EWobs.i is the measured EW, with uncertainty 
o'obs,h E^nd KWi a is from Equation[31 using the best cal- 
ibration parameters ; here the scatter parameter, (T/,q, 
is added in quadrature to the measurement error. The 
sum is over each of the three infrared passbands, a € 
{ J, H, Ks }, which have separate calibrations for the 
given spectral region (see Table [2]) . 
Since many nearby stars have 2MASS photometry, we 



E 



{g- [■M^ + 51ogio(d/10pc)]p 
2a2 



(6) 



where a corresponds to the observed infrared magnitudes 
and cTp^a is the corresponding measurement uncertainty. 

Summing over all of the indices and the additional con- 
straints gives the total of Equation [4l which we min- 
imize as a function of Mjhks, A(1/ — Ks) and d to 



6 



Pineda, Bottom, & Johnson 



determine the best-fit stellar properties and the distance 
to the star. 
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Fig. 4. — The error in the distance (top panel), color offset (mid- 
dle panel) and Mks (bottom panel) as a function of the intrinsic 
brightness of the star. We can reproduce the distances to an accu- 
racy of 11%, the color offset to 0.18 dex and Mks to 0.25 dex. The 
RMS in the absolute magnitudes for J and If band closely match 
the RMS in the K band absolute magnitude. There is no clear 
trend, as the method appears to be uniformly applicable between 
4.5 < Mks < 7.5. The solid lines mark the mean of the example 
set and the dashed lines mark the la levels about the mean. 



top panel is for A{V — Kg) and bottom for Mrs- As the 
SNR increases, the initial improvements are significant. 
But after about a SNR of ~ 70, the improvements with 
greater signal are marginal. 
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Fig. 5. — The degree to which the parameter estimates vary as a 
function of the median signal-to-noise ratio of the input spectrum 
measured across the red chip of the HIRES detector. When the 
input spectrum has a SNR above ~70 (marked by dashed line), 
the parameter estimates settle down to a well defined value. See 
text for analysis procedure. 



3.3. Accuracy 

In Figure 4 we show the results of our assessment of 
how well we can recover the properties of the stars in our 
calibration sample. The top panel of the Figure shows 
the percentage error in reproducing the observed distance 
of the stars, the middle panel shows the error on the color 
offset and the bottom panel shows the error in M-ks- 
The root-mean-scatter (RMS) for the distance, color off- 
set and absolute K-band magnitude are 11%, 0.18 dex 
and 0.25 dex respectively. Using established photomet- 
ric relations, this scatter would correspond to 0.10-0.15 
dex in [Fe/H] (dependin g on the literature calibration ) 
and ^ 0.05 in m ass (jSchlaufman fc LaughlinI [20101 ; 
iDelfosse et aLll2000t) . 

As an additional test, we examined what SNR is nec- 
essary to get consistent parameter estimates from any 
given stellar spectrum using our methods. We simulated 
a given SNR by adding noise to our template spectra us- 
ing a pseudo-random number generator. Repeating this 
many times for each SNR, we saw what effect the noise 
had on our parameter estimates. In Figure[5l we plot the 
dispersion in our estimates as a function of SNR. The 



3.4. Spectral Resolution 

We also examined how spectral resolution affects the 
utility of our calibrations. By measuring EWs consistent 
with our original measurements we would be able to use 
our calibration to recover the same set of stellar prop- 
erties using lower resolution spectra. To test this, we 
first took the template spectra for our calibration sam- 
ple and convolved the spectra down to lower resolutions 
(5,000-45,000) in increments of 5,000 from our original 
resolution of 50,000. We then normalized the spectra in 
the same manner as we did our original template spectra 
(see Section 12. 2p and computed the EW for each of the 
indices of Table [1] and all of the sample stars. In Fig- 
ure [HI we plot the average fractional difference in EW 
measurements, across all the calibration stars, between 
the convolved spectra and the unconvolved spectra as 
a function of spectral resolution, where each panel cor- 
responds to a different index as listed in Table [T] The 
error bars in the plot represent the scatter in the de- 
viation across the sample of calibration stars. For each 
index the EW measurements are consistent with the orig- 
inal measurements for spectral resolutions above 30,000. 
Therefore our calibration should not be used below this 
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threshold without accounting for the systematic effects 
demonstrated in Figure El Ahhough the integrated flux 
of a spectrum should not change at lower resolutions the 
blending of pseudo-continuum with the many absorption 
band-heads complicates our EW measurements. 
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Fig. 6. — The average fractional deviation in measured EW com- 
pared to our original HIRES spectra as a function of spectral reso- 
lution for each spectral index designated by the central wavelength 
(nm) listed in Table [l] Each point represents the mean deviation 
across the full calibration sample with the error bars given by the 
corresponding scatter. There are consierable deviations for resolu- 
tions below ~30,000. 

To consider higher resolution data we examined a sub- 
sample within our set of calibration stars that also had 
publicly available data from the ESO data archive us- 
ing the HARPS sp ectrograph at a resolution of 115,000 
(jMavor et al.l 120031 ) . For this subsample of 43 stars we 
combined multiple observations and put together tem- 
plate spectra similar to our HIRES spectra (see Sec- 
tion 12. 2|) . This produced a set of HARPS spectra with 
high SNRs, all greater than 70. We then convolved the 
spectra down to the HIRES resolution of 50,000 to com- 
pare EW measurements. Normalizing the spectra in 
the usual manner we then measured the EWs for the 
first four indices of Table [T] because of the small over- 
lap between the HARPS spectrograph and the HIRES 
red chip only these four indices were available. In Fig- 
ure [71 we compare equivalent width measurements from 
the HARPS subsample to the measurements from our 
HIRES spectral templates. The straight line in the 
plot corresponds to exact agreement. The measurements 
agreed rather well with an RMS of ^ 8%, in line with 
the combined errors between the two measurements. 

4. APPLICATIONS 

4.1. KOI-314 

We applied our methods to the case of the Kepler 
object of interest KOI-314. This M-dwarf, with a vi- 
sual magnitude 14 with a Kepler-band magnitude 
Kp — 12.93, received Keck Observatory HIRES follow-up 
to confirm the planetary nature of the two transit signals 
observed in its Kepler light curve. We were able to make 
use of 6 CPS observations to construct a high-resolution 
spectral template to apply our method. The co-addition 
procedure, as discussed in Section |2l yielded a spectrum 
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Fig. 7. — A comparison of the EW measurements for the subsam- 
ple with HARPS spectra. The different symbols correspond to the 
four different indices available and the line indicates a 1:1 corre- 
spondence. The scatter about the line is at the 8% level consistent 
with the measurement uncertainties. 

with a typical signal-to-noise ratio of ^ 250 in the red 
portion, plenty for our purposes. Figure [51 shows the re- 
sults of our analysis with the contours for each pairing 
in our five parameter fit, Mjhks, A(T^ — Ks), and d. 

Each parameter is highly correlated, leading to the 
oblong shaped contours of the figure. This is predom- 
inately because the absorption strength of the spectral 
indices is degenerate in mass and mctallicity; the ab- 
sorption strength can increase with a drop in the effective 
temperature or an increase in the metal content. Addi- 
tionally, the shape of the contours between the infrared 
magnitudes must be consistent with the observed colors 
{ J - H,H - K etc.). 
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Fig. 8. — Contour plots of the posterior probability distribution 
from the five parameter MCMC analysis of KOI-314. The contour 
levels represent the Icr, 2a, and Str confidence levels. The oblong 
contours come from spectral indices being degenerate in Mks a^nd 
/\{V-K,) 

We can then marginalize over the full posterior proba- 
bility for the five parameters to get probability distribu- 
tions for the likelihood of each of the parameters. Our 
estimates are shown in Table [31 where we have adopted 
for our uncertainties the scatter we have in reprodu cing 
the stellar parameters as demonstrated in Section 13.31 
We have additionally added an empirical estimate of the 
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TABLE 3 
Properties: KOI-314 



Attribute 


Muirhead et al. f2012a) 


This Study 


Mj 




6.18 ±0.25 


Mh 




5.57 ±0.25 


Mks 




5.39 ±0.25 






-0.20 ±0.18 


d 




66.5 ± 7.3 pc 


Mass 


0.51 ± 0.06 Mq 


0.57 ±0.05 Mq 


Radius'" 


0.48 ± 0.06 Rq 


0.54 ±0.05 Rq 


[Fe/H]<= 


-0.31 ±0.13 


-0.28 ±0.10 



^ For this work the estimate uses the relations of 
IDelfosse et al.l I I2000I) to convert Mks to mass 



For this work t he estimate uses the relations of 
IBovaiian et al.l II2012I ) to convert mass to radius 

In convertin g A(V — Ks ) to [F e/H], the estimate uses 
the relation of lNeves et "alt ll20il'l : [Fe/H] = 0.57A(y - 
Ks) — 0.17. The listed uncertainty does not include the 
scatter in their relation. 

stellar radius based on the s ingle star mass-radius rela- 
tion of IBovaiian et al] (|2012[ ) established from interfero- 
metric radii measurements of nearby low- mass stars. The 
table a lso includes select estimates from llVIuir head et al.l 
(|2012aD . which uses infrared spectra in conjunction with 
stellar models to derive properties. Additionally, in Fig- 
ure |9l we plot the marginalized distributions for the stel- 
lar properties, after converting the photometric proper- 
ties, Mks and A(y — Ks) to the physical properties of 
mass and [Fe/H]. The top left panel shows our distribu- 
tion for the distance. The top right panel s hows the dis- 
tribution for [Fe/H], using the calibration of lNeves et alj 
(|2012i) to make the conversion from A(y — Kg). The 
acc ompanying dash e d line is a Gaussian representing 
the IMuirhead et al.l (j2012af) estimate. The bottom left 
panel shows the distribution fo r mass, the sol i d line be- 
ing our estimate, employing the IDelfosse et all ()2000D re- 
lation for M.ks with t he dashed line representing the 
HVluirhead et al l (|2012aD estimate. Finally, the bottom 
right pane l uses our mass estimate and the mass-radius 
relation of IBovaiian et al.l (j2012D to estimate the radius 
of the star with t he das h-dot line again representing the 
HVTuirhead et al.l (l2012af) estimat e. Our measurements 
match those in llUuirhead et al.l (2012a) within the re- 
spective uncertainties, giving us confidence in the accu- 
racy of our derived parameters. Combining our estimates 
with theirs gives a mass estimate of 0.55 ± 0.04 Mq, a 
radius estimate of 0.52 ± 0.04 Rq and a metallicity of 
-0.29 ±0.08 dex . 

4.2. G J 3470 

As an additional application of our methods we con- 
sider t he exoplanet hosting star GJ 3470. In the discovery 
paper iBonfils et all ()2012( ) used HARPS radial velocities 
and a photometric transit detection to constrain the mass 
and radius of the planet, adding the system to three other 
M dwarfs (GL436, GJ 1214 and KOI-254) with planets 
that have well-measured mass and radius estimates. The 
precision of the planet properties for GJ 3470b however 
was limited by the uncertainty in the stellar properties. 
Previous studies of GJ 3470 suggested that it is a typical 
field star on the main sequence, making the methods of 
this paper applicable ()Bonfils et al.l[2012l ). 

In Section 13.41 we showed how HARPS archival spec- 
tra could be used to measure EWs compatible with 
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Fig. 9. — Marginalized probability distributions for the proper- 
ties of KOI=314, using photometric calibrations to convert to mass 
and [Fe/ H]. TheQP clfossc ct al. (2000) is used to convert to mass 
and the INeves et a l. (2012) relation is used for metallicity. The 
distributions from Muirhead ct al. (2012a) are over plotted with a 
dashed line, assuming normal distributions with a standard devia- 
tion given by the reported uncertainties 



TABLE 4 

Properties: G.J 3470 System 



Attribute 


Bonfils et al. ('2012) 


This Study 


J'' 

a/iJ* 
Period 
Rp/ R-k 


14.9 ± 1.2 
3.33714 days 
0.0755 ± 0.0031 


8.794 ±0.019 
8.206 ± 0.023 
7.989 ± 0.023 


d 

Mass 
Radius 
[Fe/H] 


0.54 ± 0.07 Mq 
0.50 ±0.06 Rq 


29.9±l l pc 
0.53 ± 0.05 Mq 
0.50 ±0.05 Rq 

0.12 ±0.12 


2MASS Dhotometrv oflCutri et al.l (|2003ll 



the spectral index calibrations of Table [5] We applied 
the procedures of Section 13.41 on the HARPS spectra of 
GJ3470 to obtain a spectrum with a typical signal-to- 
noise of 72 near the spectral indices. We measured the 
EWs and were able to independently estimate the dis- 
tance, mass, [Fe/H] and radius of the star using entirely 
empirical methods: w e used the mass-radius relation of 
iBovajian et al.l (|2012( ) to determine the radius from the 
mass. However, only the first four of the indices of Table 
[T]are useable with the HARPS spectra, reducing the pre- 
cision of our measurements. Using just the four indices 
we reproduced the properties of the calibration sample 
to a RMS of 0.38 in Mjhks, 0.27 dex in A{V-Ks) and 
18% in distance (see methods in Section 13731) . This RIVIS 
corresponds to 0.08 AIq in m ass and 0.15 de x in m etal- 
licity using t he cal ibrations of IDelfosse et al.l ()2000[ ) and 
INeves et al] ()2012D respectively. To improve the preci- 
sion of our estimates we look for additional constraints 
to apply to the stellar properties. 

The broadband photometric methods of lJohnson et al.l 
lf2012a) showed how photometric observations could be 
combined with transit light curve observables to provide 
precise estimates of the stellar properties. Following their 
example we include the reduced semi-major axis of the 
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planet orbit as an additional constraint: 



(-) 

U7r2 J 



1/3 



1/3 



_p2/3 



(7) 



where the radius is giv en as a function of ma ss using the 
mass-radius relation of lBovaiian et all ()2012D and we ne- 
glect the mass of the planet as very small compared to the 
stellar mass. The transit light-curve observable, dp/i?^, 
gives the planet-star separation at the time of transit 
and for the case of a circular orbit matches the reduced 
semi-major axis of Equation [T] We take the period to be 
well defined as 3.33714 days and the observed (a/i?^)o 
as 14.9 ± 1.2 from Bonfils et al. (2012). We incorporated 
this constraint by adding an additional term to the total 
of Equation [J 



(8) 



In implementing t he MCMC methods of Section 13.21 we 
use the relations of lDelfosse et al.l (j200ClD to go from Mks 
to M-i, in Equation |51 Our estimates, shown in Table |31 
agree well with the output values derived by the planet- 
discovery team which combined a transit detection with 
radial velocity measurements, 0.54 ± 0.07 Mq for the 
mass and 0.50 ± 0.06 Rq for the radius, however our 
stellar properties are more precise. They are also con- 
sistent with stellar parameter estimates incorporat ing a 
new infrared transit analysis (jPemorv et al.l 12013 '). In 
Figure [10] we plot the contours for the total probabil- 
ity distribution in our estimate of the stellar mass as a 
function of observed quantities at the fixed distance and 
metallicity of our estimates shown in Table |31 The filled 
overlays show the 3ct constraints provided by each ob- 
servation with blue corresponding to a/i?*, green to the 
joint constraint in JHK and orange given by the com- 
bined constraints of the EWs. 

We can combine our stellar radius m easurement with 
the estimate from IBonfils et al.l (|2012D to get a precise 
stellar radius of 0.50±0.04 Rq. Usi ng the transit obsery - 
able Rp/R^ = 0.0755 ± 0.0031 from ' Bonfils eFall (IMl) . 
we get a planet radius estimate of Rp = 4.12 ± 0.37 i?®, 
slightly smaller than their estimate of i?p = 4.2 ±0.6 i?®. 

5. SUMMARY AND DISCUSSION 

Using high-resolution spectra of nearby M-dwarfs we 
have developed a new spectroscopic calibration for the 
physical properties of low-mass stars of types K7-M4. 
With a high signal-to-noise spectrum and 2MASS pho- 
tometry we can estimate the mass and the metallicity of 
the star as well as the distance to the star. 

Our methods are based on the integrated spectral flux 
(EW) at a series of spectral regions sensitive to the 
known photometric properties of the stars. By mea- 
suring the corresponding EWs from the spectra of low- 
mass stars and making use of the observed infrared mag- 
nitudes, we can estimate the photometric properties, 
■MjHKs and A{V — Kg), as well as the distance. Ad- 
ditionally, using known photometric calibrations, these 
can be converted to mass and metallicity. Our estimates 
are strongly dependent on the accuracy of the assumed 
parameters for our calibration stars, however the sample 
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Fig. 10. — Contours of the probability distribution for our mass 
estimate of GJ 3470 as a function of observed parameters with the 
metallicity and distance fixed to their best estimates, 0.12 dex and 
29.9 pc respectively. The outer (blue) region defines the Scr con- 
straint applied by the measurement of a/_R*. The middle (green) 
regions applies the combined 3cr constraint of the infrared pass- 
bands JHK. The thin (orange) region is the combined 3(t con- 
straint using our calibration on the measured EWs of the HARPS 
spectra. The black contours represent the Icr, 2(t, and Scr contour 
levels for the entire set of constraints. 
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Fig. 11. — Contours indicate the approximate minimum total in- 
tegration time with HIRES, in minutes, for typical low-mass star 
required to utilize our methods to obtain stellar properties (achieve 
S/N ~ 70). The time is shown in the black labels and the corre- 
sponding approximate V-hand magnitude is shown with the red 
labels 

comprises a set of well studied nearby stars with accu- 
rate properties. We are thus able to estimate the intrin- 
sic physical properties of low-mass stars without having 
parallax measurements, independently of stellar models. 

We applied our methods to a particular star KOI-314, 
one of the Kepler objects of interest and estimated a 
mass of 0.58 ± 0.05 Mq, a radius of 0.55±g:[;| Rq, a 
metallicity, [Fe/H], of —0.28 ± 0.10 and a distance of 
66.5 ± 7.3 pc, where we have adopted for our uncertain- 
ties the representative scatter of Section 13.31 and propa- 
gated those uncertainties through the empirical calibra- 
tions to the uncertainties of the desired physical proper- 
ties. These estim ates are in good agreem ent with the ap- 
proach taken by iMuirhead et al.l ()2012aD , providing ad- 
ditional observational evidence in corroboration of their 
model-dependent methods. We were also able to ap- 
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ply our methods to the planet host G J 3470, making use 
of archival HARPS spectra and taking advantage of the 
transit observable a/ Ri, to narrow in on the stellar prop- 
erties. We estimated a mass of 0.53 ± 0.05 M©, a radius 
of0.50±0.05 i?0, ametahicity, [Fe/H], of 0.12±0.12 and 
a distance of 29.9±3 4 pc. These propert i es are very simi- 
lar to those determined bv lBonfils et al.l (|2012[ ) who then 
estimate an approximate mass of '^14 and radius of 
~4.2 i?e for GJ 3470b. 

Although our calibration provides estimates for stellar 
properties, there are some limitations. We examined the 
effect of spectral resolution and found that we could use 
the calibration only for spectral resolutions greater than 
'^30,000. We also demonstrated how to use the calibra- 
tion with higher resolution data (>50,000) by making 
use of HARPS archival spectra. Additionally, our cal- 
ibration sample only spans a particular range of Mks 
(see Figure[T]) and should only be used outside that range 
with caution. This restricts applicability to mostly early 
type M dwarfs, earlier than about M4, and late K dwarfs. 
This still spans a fairly broad range of masses from about 
- 0.7 to - 0.2 Mq. 

Our method requires a high-signal, high-resolution, 
spectrum of the star. To quantify the necessary signal- 
to- noise, we used the observations of KOT314 as a guide 
in our noise analysis (see Section [373]) . In Figure [TTl we 
show the approximate minimum total integration time 
needed, in minutes (on left side), to achieve a signal- 
to-noise of ^ 70 and use the techniques presented in this 
contribution. The corresponding V-band magnitudes are 
also shown on the contours towards the right side of the 
plot. As a benchmark, a star with a l^-band magnitude 
of 14 would need a total integration time of 30 minutes. 
It is possible to build up this signal over time by building 
a composite spectrum. This makes it an ideal method to 
complement radial velocity surveys of M dwarfs. Many 
spectra are needed to sample the radial velocities of these 
stars, so as a byproduct of those observations the phys- 
ical stellar properties can be determined simultaneously. 
Additionally, this method can also be applied immedi- 
ately to archival HIRES data of low-mass starfl 

APPENDIX 
CONTINUUM NORMALIZATION 

The continuum normalization procedure was briefly 
explained in Section 12.21 however, we expand upon the 
details here. The reduced spectra from the HIRES de- 
tector include several orders, ten of those corresponding 
to the red chip of the detector, spanning between 650 
nm and 800 nm. Each spectral order is affected by the 
blaze function of the detector and the overall shape of the 
stellar spectrum, and we therefore normalize each order 
separately. For each order, the spectral regions listed in 
Table [SJ we first masked out any telluric regions. We 
then equally divided each region into 10 or 20 (see 'Divi- 
sions' in Table [5]) different sections. For each section we 
ordered the flux and took the top 1% or 2% (see 'Per- 
centage' in Table [5]) level as representative of a pseudo- 
continuum. We then fit these points with a polynomial 
of order 2 or 3 (see 'Polynomial Order' in Table [S]) to 
define our normalization. The number of points in the 

^ Ke ck Observatory Archive: 
|http://www 2. keck. hawaii.edu/ koa/public/koa.php I 
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Fig. 12. — Example spectrum with normalization plotted over- 
top. Each plot comprises a full order and shows the regions used as 
spectral indices for computing equivalent widths. Points used for 
the continuum normalization are shown as filled circles. Chosen 
points assure that the normalization addresses the convex shape 
of instrument profile. All of the wiggles are real features in the 
spectra. 

fit therefore matches the number of bins we used. For 
each order we divided the spectrum by this polynomial 
to get the normalized spectrum (see example in Figure 
I12p . Changes in the normalization properties attempt 
to account for differences in the curvature/symmetry of 
the spectral continuum/blaze function to achieve an ap- 
propriate normalization. To this end the fifth spectral 
order in Table [5] does not use the full order in the nor- 
malization and instead we applied a narrower range of 
wavelengths when applying our normalization methods 
in order to account for the broad and deep absorption 
bands (see bottom panel of Figure [T2|). Despite the lack 
of pseudo-continuum points toward greater wavelengths 
in the lower panel, the chosen points assure that our nor- 
malization accounts for the generally convex shape of the 
blaze function for each the order. Although there can be 
issues at the edge of each order, the spectral indices are 
all located near the central regions so this does not affect 
the calculated equivalent widths. 

We opted not to use continuum regions defined as the 
linear interpolation of points flanking the absorption re- 
gion because the spectral regions we identified as sen- 
sitive to stellar properties were not always bounded by 
points of minimal absorption. Our procedure also ad- 
dresses the difficulties introduced by the shape of the 
blaze function on the spectrum and by taking a broader 
region into account for the continuum calculation we 
can get consistent pseudo-continua across many differ- 
ent spectra. Our application of this procedure to the 
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TABLE 5 
Normalization Properties 



R fieri nn ( ri m 1 

±Lt^tiIl_Jll ^ 11111 J 


1_-/1 V lolWllo 


X V^l ^-'V^llLditiL^ /{) 


Mnl ^mr\m 1 a 1 (ivrifiv 
1 uiyiiuiiiidii v_y 1 i_n.,i 


654.35 - 665.92 


10 


1 


3 


666.69 - 678.49 


10 


1 


2 


679.52 - 691.53 


10 


1 


3 


692.85 - 705.08 


10 


1 


3 


706.71 - 719.17 


10 


1 


3 


721.14 - 733.84 


10 


1 


3 


736.17 - 749.12 


10 


1 


3 


751.84 - 765.05 


20 


2 


2 


768.19 - 781.67 


10 


1 


3 


785.28 - 796.55 


10 


1 


2 



^ Within region only the following sections are used for computing 
normalization: 707.04-708.79 nm, 712.19-712.58 nm, 718.35-718.63 



HARPS data was able to produce equivalent width mea- 
surements for several indices in accord with our HIRES 
measurements, showing how our methods can be used for 
different data samples (see Section [3^ . 
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